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Abstract 

We numerically perform a spectral analysis of a quasi-periodically driven 
spin 1/2 system, the spectrum of which is Singular Continuous. We compute 
fractal dimensions of spectral measures and discuss their connections with 
the time behaviour of various dynamical quantities, such as the moments 
of the distribution of the wave packet. Our data suggest a close similarity 
between the information dimension of the spectrum and the exponent ruling 
the algebraic growth of the 'entropic width' of wavepackets. 

I. INTRODUCTION 

The increasingly frequent apparition of Singular Continuous (S.C.) spectra in various 
quantum mechanical situations has attracted attention to the dynamical implications of such 
spectra. For example, for an electron moving in an incommensurate or disordered structure, 
singular continuous spectra typically result in a sort of pseudo-diffusive dynamics, which 
has a direct bearing on transport properties. If the spectrum is a multifractal, some of these 
properties depend on the value of certain fractal dimensions, according to quantitative rules, 
the precise assessment of which is an important theoretical task . 

SC spectra have also been found in some periodically or quasi-periodically driven model 
systems, which were introduced in order to investigate the possibility of chaotic behaviour 
in quantum dynamics. [@] @. In particular, quasi-periodically driven spin systems have 



been studied, in view of their formal simplicity. Various spectral types have been identified 
, on varying levels of evidence, for different systems in this class 0] . From the mathematical 
viewpoint, some quasi-periodically driven spin systems formally belong to a class of abstract 
dynamical systems, which is well-known in ergodic theory, and for some of which SC spectra 
have been rigorously proven to occur |^. In a recent paper |^ a renormalization group 
analysis has been implemented, strongly supporting SC spectra for a wider class of spin 
systems than encompassed by available exact results. The formal simplicity of this class of 
systems makes a numerical analysis of their dynamics particularly convenient, so that they 
appear very well suited for the study of the dynamical implications of fractal spectra. 

In this paper we describe a numerical spectral analysis for a particular model in this 
class, aimed at computing certain fractal dimensions, and at connecting them to asymptotic 
aspects of the dynamics. We use a technical approach, based on a discrete-time variant of 
Floquet theory, which has not been implemented before in this context, and which offers 
a twofold advantage. In the first place, if combined with a suitable scheme of rational ap- 
proximation for an incommensuration parameter, it allows for a reliable computation of the 
spectral measures, at a relatively low computational cost. Second, the Floquet formulation 
allows for a meaningful analysis of the growth in time of the spread of wave-packets, which 
is a very efficient empirical marker for continuous spectra, and usually a more convenient 
one than the hitherto used decay of correlations. 

We then analyze a number of dynamical features, related to the growth of wave-packets. 
We find that the exponents of growth of the moments of the wave packet are not directly 
related to the Hausdorff dimension of the spectrum, (which is 1, as shown by a simple 
rigorous argument). We also find that the exponent of growth of the "entropic" spread of 
the wavepacket is close to the numerically computed information dimension of the spectral 
measure. Finally, we shortly present a general bound for the growth of the information 
contents of a finite string of observations with the length of the string, in terms of the 
information dimension of the spectrum. 
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II. THE MODEL 



Our model is a periodically kicked spin | system, with the kicking strength depend- 
ing quasi-periodically on time. We consider a function S{ip) from [0, 2tt] into the unitary, 
unimodular 2x2 matrices, given by: 

S{^) = e^'=^('^)*^e*^^ (2.1) 

where x is a periodic real- valued function, to be specified later, /c is a parameter, and ax, Oz 
are Pauli matrices. In [0, 27r] we consider the shift Tq : </? i — > </? -|- 27ra{mod27r) , with a a 
fixed parameter. We define the dynamics of our model by first arbitrarily fixing a phase (po, 



and then prescribing that spinors = 
t + 1 into: 



given at (discrete) time t evolve at time 



V'(i+l) = ^(T>o)#) (2.2) 

At fixed ifio, a, k the evolution of a given initial spinor is thus obtained by applying a sequence 
of unitary matrices, which can be either a periodic or a quasi-periodic one, depending on 
the arithmetic nature of the number a :if the latter is a rational, a = p/q with q mutually 
prime integers, the sequence is periodic with period q, otherwise it is quasi-periodic. We 
have chosen the function x ^ the periodicized characteristic function of an interval / in 
[0, 27r] : xi^) = 1 if (p{mod2-K) E I, %((/?) = otherwise. As we shall explain below, such a 
choice makes the numerical analysis very efficient. The sequence of unitary matrices defining 
the evolution is now uniquely defined by the symbolic trajectory of the chosen (po associated 
with the shift t„ and with the partition of [0, 27r] defined by / and its complement. Since this 
partition consists of two sets, the symbolic sequence can be written as a binary sequence. In 
particular, if the length of the interval / is taken 2Tra and a , and a is the (inverse) Golden 
Ratio (^(5) — l)/2, one obtains the Fibonacci sequence. 



III. SPECTRAL MEASURES. 



There are two possible approaches to the spectral analysis of the system ( |2.1| ). The first 
consists in considering (O) as a classical dynamical system, the state of which is defined 



by a couple 0, ip where is an angle and ip is a. normalized spinor; the one-step evolution of 
the system is given by: 

— * ^ — * 

(y?, ip) I — ^ {Ta^,S{^)i)) (3.1) 

The system thus defined belongs to a class of dynamical systems known as "skew- 
products" Spectral measures are defined from the Fourier transform of correlation 
functions; in particular, given an initial spinor "0(0), the corresponding spectral measure d^x 
is defined by: 

it\. 



R(t) = I e^*^d/i(A) (3.2) 
where R{t) is the correlation function: 







R{t) =lim lf](^>),.0(s + t)) (3.3) 

where (, ) is the scalar product in C^. Spectral measures defined in this way will be termed 
"dynamical" measures in the following. 

Another definition of spectral measure rests on a generalization of Floquet theory, which 
is obtained on imbedding the nonautonomous quantum dynamics |2.2| in an autonomous 
dynamics, defined in a larger Hilbert space. This is done as follows. We consider the phase 
as a new dynamical variable, and thereby consider state vectors \E' =ip{(p) in the Hilbert 
space H = L'^{[0, 2n]) (g) C^. In this space we consider the discrete unitary group generated 
by the unitary Floquet operator: 

(5V^)(v.)=5(^)^(r-V) (3.4) 

If '?/'(v?) is regarded as a curve in the phase space of the dynamical system (|3.1j), then 
( p.4|) specifies the evolution of this curve under the dynamics ( |3.1|) . Our second definition 



of a spectral measure is just the usual one for unitary operators, applied to the operator 
( p.4| ). Such spectral measures will be termed "Floquet measures". The connection between 
dynamical and Floquet measures will be discussed separately for the incommensurate and 
the commensurate cases. For the time being, we shall point out a simple property of the 
spectrum of S, which is valid in both cases. On defining in Ti. the unitary operator Uip{if) = 
e^'^ip{ip) we find: 

U^SU = e-2™5 

which implies that the spectrum of S is invariant under the shift : that is, if a point of 
the unit circle belongs in the spectrum of S, so does the whole orbit of that point under . 

It is worth remarking that the operator S can be interpreted as the Floquet operator for 
a linear kicked rotator endowed with spin; the occurrence of SC spectra for a spinless linear 
kicked rotator has been discussed in [§. 

IV. THE COMMENSURATE CASE. 

In the commensurate case a = p/q, the model is a periodically driven one, and its 
dynamics can be understood, as usual, from a spectral analysis of the one-period evolution 
operator. The state vector at times multiple of the period q is found from: 

^finq) = f"{^o)m 

where 

T(^o) = n^(r>o) (4.1) 

t=0 

(the product being ordered from right to left) yields the evolution over one period. This 
operator does not depend on time; thus the dynamics, read at integer multiples of the period, 
is given by a discrete unitary group, generated by a fixed unitary (|4.1| ). Being a finite matrix, 
the latter has obviously a pure discrete spectrum, consisting of two conjugate eigenvalues 
Zi = e*'^, Z2 = zl = e"*'*'. The dynamics is therefore recurrent. 
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From QdJ) we find: 

iSy){cp)=f{if)^f{ip) (4.2) 

that is, iS'' is a fibered operator. Its fibers are precisely the 2x2 matricesT((y9), and its 
spectrum is the range of the functions zi{if ),Z2{(f ) giving the eigenvalues of T{if ). It turns 
out that, when the size of the interval I is not a multiple of 27r/g, both functions have 
exactly two values in their range: otherwise, they have just one. This fact allows for a 
decisive simplification of the numerical analysis, and is proved in Appendix A; it stems from 
the finiteness of the range of the function x- 

Therefore the spectrum of iS^ consists of a finite number 2M of eigenvalues zjsf = 
e'^^^,{N = 1,...2M), in complex conjugate pairs, with M = 2 or M = 1. The spectrum 
of S is then a subset of the set of the 2Mq complex g-th roots of the z^^. The corresponding 
eigenphases are: 

Am 2TTj , , 

Xj,N = ^ + — ^ 4.3 

The spectrum of S is discrete, i.e., the Floquet dynamics in the extended Hilbert space has 
the same spectral character as the proper dynamics. It should be noted that the fact that 
the dynamical measure and the Floquet one are of the same type is, in the commensurate 
case, a nongeneric feature , connected with the choice of a finitely valued function x in ( P-2[) . 
On choosing a nonconstant analytic function x, an absolutely continuous spectrum of S 
would be found instead, in contrast to the character of the dynamics, that would be still 
recurrent. 

V. THE INCOMMENSURATE CASE. 

The spectral measure of a vector =ip{^p) with respect to the operator ( |3.4| ) is defined 
as the Fourier transform of the correlation: 



C{t) = {^,S'^)^ = / rf¥.(V^M,^M^(r-V)...5(r>Mr-i>)) (5.1) 
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In the incommensurate case, the well-known ergodic property holds: 

If 1 ^ _ 

— / dy^g{y^) = lim - ^ g{T^ Vo) (5.2) 
27r J T^oo T 

for any summable function 5f,and for almost every (fo. Then the integral over (p in ( |5.1| ) can 
be computed as follows: 

C{s) = ^lim i g (^(cVo), ^(r~Vo)^(r„'*'Vo)...^(r„-*-Vo)^(r-^-*-Vo)) (5.3) 

Now we apply to both factors in the scalar product the unitary operator 
S{ipo)S{T~^Lpo)...S{T^^^^ipo) and take ip = ip{0). Then, on comparing the result with ( p.3|) 
we find that 

C{t) = R{t) except possibly for a set of zero measure of values of ipQ. Therefore, in 
the incommensurate case the dynamical spectral measure ( |3.3| ) coincides with the Floquet 
spectral measure. 

In a quite similar quasi-periodically kicked spin system was studied, which was reduced 
to a classical dynamical system for which singular continuity of the spectrum is a known 
mathematical result. It looks likely, on account of the closeness of our model to that one, 
that this result holds in our case, too; our numerical results, presented below, yield strong 
evidence in this sense. 

Since in the incommensurate case the orbit of any point of the unit circle under the 
shift is dense in the circle, from the invariance of the spectrum under Tq, it follows that the 
spectrum (which is by definition a closed set) must coincide with the whole unit circle. 

VI. NUMERICAL RESULTS 

Our main focus is on the incommensurate case, with a given by the Golden Ratio. In 
most of our computations we have taken / an interval of length 27r/3; in a few cases we 
have taken 27ia instead. Unless explicitly stated, we will make reference to the first choice. 
We have obtained two types of numerical data: dynamical ones, obtained from directly 
simulating the Floquet dynamics, and spectral ones. We discuss the latter first. 
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A direct computation of the spectral measure from the Fourier transform of the nu- 
merically computed correlation function ( |3.3| ) involves certain subtleties, well-known in the 
spectral analysis of time series, connected with the necessity of appropriately weighting 
the tails of correlation functions [|I0|. A more convenient strategy for the present model is 
based on approximations of the Golden Ratio constructed via continued fraction expansion, 
which produces a well-known sequence of rational approximants. For such approximants we 
carefully compute the spectral measures, as we explain in Appendix B, and thus we obtain 
approximations of the true spectral measure, by means of the pure point measure ( |9.2|) . 

A similar strategy has been widely used in the numerical investigation of other incom- 
mensurate spectral problems, such as e.g. the Harper and the Kicked-Harper model ||^. It 
is worth remarking, however, that in those cases, unlike the present one, the commensurate 
approximations have an absolutely continuous band spectrum, that is usually taken as a 
covering of the limit spectrum. 

For the obtained measures we have performed a fractal analysis, aimed at the determina- 
tion of the fractal dimensions Dg] the Hausdorff dimension Dh is 1 in the incommensurate 
case, because the spectrum in that case is the whole unit circle. In our computations we 
have used the definitions, 

Di = lim — ; ; L'„ = lim- — — - (6.1) 

<5-o In 5 ^ .5-0 (g - 1) In 5 ^ ^ 

where the interval [0, 27r] was partitioned into small intervals of size 6, the i—th of which 
receives a weight fii from the spectral measure. 

Dq were obtained by a linear fit of the numerators in ( |6.1|) vs the denominators in a 
suitable range of 6. A typical result is shown in fig.l. In no case could the value -Dq obtained 
from the spectrum of generalized dimensions be distinguished from Dh = 1. The correlation 
dimension D2 is known to rule the decay of integrated correlations 1^: 

^ElCm'- const.. T-^^ (6.2) 

asymptotically for T — > 00. We have used this exact result to adjust our numerical com- 
putation of the other dimensions; in fact, since the approximating measure is a pure-point 
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one, the dimensions of the hmit measure must be read off a suitable range of not too small 
6. Our appreciation of a reliable range for 6 was based on the comparison of the rhs of the 
2nd eq.( |6.1D with the exponent 7 of algebraic decay of correlation obtained from a direct 
computation of the latter (Fig. 2). 

We have also computed a number of dynamical data, by a direct numerical simulation 
of the dynamics. In particular, we have analyzed the growth in time of the momenta la of 
the distribution of the Floquet wave-function on the Fourier basis. /„ are defined by: 



s=Oj=l 







(6.3) 



Notice that moments of order a > 1 diverge for t > 0, because the Fourier coefficients of 
the wave function behave as at large n, due to the discontinuity of x ■ 
Finally we have studied the growth in time of the average entropy S{t) : 

+00 

n=— 00 

Examples of the dependence on time of the moments and of the entropy are given in 
figs. 3,4. The moments increase according to a power law: /^(t) ~ const. t^^'^°'\ and so does 
the "entropic number of states": N{t) = exp(5'(t) ~ const.f^ . In Table 1 we summarize the 
values of the exponents 7,/5, a and the values of the Di, D2 for several choices of the kicking 
strength k. 

The specified errors are those involved in linear fits of bilogarithmic plots used to compute 
fractal dimensions or growth exponents; they do not include numerical errors in computing 
the data themselves, therefore they to some extent underestimate the real errors . Generally 
speaking, it is usually difficult to get precise estimates for the dynamical growth exponents 
associated with fractal spectra, partly because of finite-basis effects coming into play at 
large times, and much more because the curves of growth display a characteristic pattern of 
kinks . In cases in which such seemingly log-periodic structures were particularly evident, 
the growth exponent was obtained directly by drawing a straight line through a sequence of 
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maxima. In such cases the fitting error, not specified in Table 1, is on the order of the first 
missing digit. 



VII. DISCUSSION. 

Some general exact results are known about the long-time properties of the quantum 
dynamics in the presence of a SC spectrum. First of all, integrated correlations like ( |6.2|) 
must tend to zero, and momenta must diverge in the limit t ^ oo in all cases when the 
spectrum is purely continuous; furthermore, if the spectral measure is a fractal one, with a 
correlation dimension D21 and information dimension Di, then (|6.2| ) holds for the decay of 
correlations 0, and (3{a) > -Di 0. Data in Table 1 are fully consistent with these exact 
estimates (as mentioned above, the estimate ( |6.2| ) was actually used to adjust our numerical 
method). In addition, more or less heuristic arguments relating more precisely the exponents 
(3{a) to multifractality have been attempted. One such argument ||^ has led to hypothesize 
that (3{q) = Dh, the Hausdorff dimension of the spectrum. Although consistent with 



numerical results from a number of models, this hypothesis has been called in question [jTT 
and also recent numerical investigations have provided evidence that (3{a) is not in general 
a constant but covers a continuous range of scaling exponents (multi-scaling). Our present 
data are fully consistent with the latter picture, because the observed values of (3{a) differ 
from Dh = 1 (which was theoretically established and numerically confirmed) significantly 
more than the estimated numerical error. Somewhat smaller, but still significant, are the 
differences observed in the values of /3(a) obtained for different values of a. 

Two interesting facts emerging from table 1 are (i) the closeness of the values obtained 
for different values of k, which seems to indicate that the fractal structure of the spectrum 
is essentially determined by the quasi-periodic structure of the symbolic strings alone, and 
(ii) the closeness of the value of the "entropic" exponent a to that of Di. Although the 
difference of the two values was comparable to the relatively large numerical error, the 
actual agreement may be much better, because many of the data used in estimating Di 
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appear not to have fully converged to the proper asymptotic regime. These data are in fact 
compared to the dynamical entropic exponent in Fig. 5 ; the agreement "by eye" is there 
better than implied by Table 1. Whether this fact reflects a real connection between the two 
quantities is an interesting theoretical question, for which we have no answer for the time 
being. 

Finally we shall mention a dynamical property, that, although not directly related to 
the above numerical results, has been recently introduced in essentially the same context as 
discussed here, with the aim of analyzing the possibly chaotic property of the evolution of 
the spin system. For a classical dynamical system in discrete time, with a compact phase 
space Q, the "informational complexity" is introduced by considering a finite partition 11 
consisting of m subsets, and the associated symbolic dynamics. For any given integer time 
t, all the possible symbolic strings of length t define a partition 11* of If a specific orbit 
of the system is chosen, then a probability can be attached to every class of this partition, 
defined as the frequency with which the given finite string appears in the infinite symbolic 
string of the given trajectory. The Shannon entropy H{t) can then be defined, as well as 
the corresponding number of histories N{t) = exp{H{t)). A chaotic behaviour is associated 
with an exponential growth of N{t) in time. 

In a recent paper [l^ such an analysis has been numerically implemented for the case of 
a quasi- periodically driven spin system. In essence, Q was taken as the compact variety of 
normalized spinors; a certain partition of Q was constructed, and the growth of H{t) with 
time was analyzed. A seemingly sub-exponential growth of N(t) was observed, of the type 
exp(ct'^) with 7 < 1. 

We shall here sketch a general argument, which establishes an upper bound for N{t) in the 
presence of a fractal spectrum. At time t, 11* is a partition of a t. //-dimensional space, where 
u is the dimension of Q; therefore H{t) < tulogm. Now, all the strings iIj{s + 1), ...ipls + t) 
of length t which are observed in the evolution do indeed span a vector space of dimension 
t, because of the pure continuity of the spectrum; however, their "effective dimension" can 



be much less. It can in fact be proven |T3| that the minimum dimension d^{t) of a subspace 
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which contains all the strings of length t within a maximum error e asymptotically grows 
like This means that, apart from a small error of order e, the entropy H{t) is the same 
as the one computed over ~ m^'^^^^ classes. This leads to the upper bound H(t) < const.t^^. 

Generally speaking, a direct numerical analysis of N{t) seems difficult, because comput- 
ing the frequency of strings long enough to reproduce the correct asymptotic regime requires 
an enormous computation time. In any case, for the special case of quasi-periodically driven 



spin systems, the above upper bound is but a very crude one: in fact a simple argument ||I4 
indicates that the asymptotic growth of N{t) cannot be faster than algebraic. 



VIII. APPENDIX A 

The matrix T{(f>) is uniquely associated with a binary string of q digits, which specifies 
the symbolic periodic trajectory of if. Suppose that the strings corresponding to two matrices 
T{ip),T{if') differ from each other merely by a cyclic permutation of digits. Then the two 
matrices differ from each other by a cyclic permutation of the S operators which enter as 
factors in the definition of the matrices themselves. The two matrices are then unitarily 
equivalent, and have therefore the same eigenvalues. We have thus reached the conclusion 
that the number of distinct values in the range of Zi{ip) (and of -22 (v^), as well) is equal 
to the number of nonequivalent symbolic strings of q digits, two strings being equivalent if 
they can be obtained from each other by a cyclic permutation. Thus in order to find this 
number we have to find the total number of symbolic strings of length g, and to divide it 
by g, which is the number of strings in an equivalence class. Let us denote by Co and Ci 
(with Ci = I) the two classes of the partition. Then the points (f which produce a given 
symbolic string ii,i2, ■■■iq are those belonging to the set Ai-^i^,,,i^ = r^^^Ci-,^) fl ... fl T~^(Cjg). 
The sets Ai-^i^,,,i^ define a partition of [0, 2tt] , and there are as many distinct symbolic strings 
of length q as are nonempty classes in this partition. It is easily seen that there are either 
2q or q nonempty such classes, the latter case occurring if, and only if, / is a multiple of 
27r/g; in fact, on repeatedly applying the shift Tq, to the interval J, we get q distinct intervals, 
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whose endpoints make a set of 2q or q distinct points, depending on whether the length of 
J is a multiple of 27r/g,or not. The sets Aj^jj. are precisely the disjoint intervals in which 
[0,27r] is divided by these 2q (resp., q) points, and their number is therefore 2q (resp., q). 
The number of nonequivalent strings is thus exactly 2 (resp., 1). 



IX. APPENDIX B 



In a commensurate case with a = p/q we get, from the definition :: 



VI/) = ^ e 



2M q-1 



N=l 



j=0 



^2mtjq 1 



H 



(9.1) 



where Pj^n is projection over the eigenspace of c corresponding to the eigenphase \j,N- The 
spectral measure in [0, 2tt] is pure point: 



2M g-1 2 

EE||4iv^L'5(A-A,,;v) 



(9.2) 



N=l j=0 

In order to compute this measure we have to find the eigenphases Xj^n and the associated 



weights pj^ 



N 



H 



. The eigenphases are most easily computed, by diagonalizing the 



2x2 matrix (|4.1|) , for two different values of ipo (e.g., the two extremes of the interval /), 
and then using ( |4.3| )(we assume that q is not a multiple of 3, for in that case just one value 
of (po would be sufficient.). Concerning the weights, denoting P/v = J2jf'j,N, we have from 



-itAisrq 



q-1 



(9.3) 



j=0 



which is an exactly periodic function of the discrete time t with period q, so that the weights 
are quite easily and reliably computable via finite Fourier transform as soon as the Ihs of 
( |9.3D is known. Since A^r is computed as described above, we are left with the computation 
of the correlation of Pn"^ at times t = 0, ...q — 1. To this end we first find -FW^ as follows: 
having diagonalized the matrix ( [4.1D for all values of cpo in a suitably thick grid, we let 
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where UN{fo) is defined as follows: it is the eigenvector of ([4.1|) with eigenvalue A^, if 



the latter is an eigenvalue of ( |4.1| ) at the given (po, and it is the zero vector otherwise. 
The correlation of Pn"^ is then found directly from the definitions (|3.4| ),( pTTD . The only 
approximation involved in this computation is the discretization of the scalar product in Ti, 
which involves an integral over (p and is instead computed as a finite sum over the chosen 
finite grid used to discretize [0, 27r] . 



We thank R.Artuso for valuable advise in computing fractal dimensions, and A.Vulpiani 
for an useful discussion. 
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X. FIGURE CAPTIONS 

Fig.l - Spectrum of generalized fractal dimensions Dq, for k = 10.5. The rational ap- 
proximant of the Golden Ratio used in this computation was 1597/2584. 

Fig. 2 - A bilogarithmic plot of J2 yU? versus (eqn. (|6.ip ?he slope of the straight line is 
the exponent of the correlation decay, eqn. (|6.2|) . Here k = 10.5. 

Fig. 3 - Illustrating the growth of the moment I1/2, for k = 4. 

Fig. 4 - Illustrating the growth in time of the entropy for k = 10.5. The slope of the 
straight line is the computed information dimension Di. 

Fig. 5 - A bilogarithmic plot of — Z^yUjlog/Xj versus (eqn.(|6.1|), for the case = 1 of 
Table I. The slope of the dashed line is the entropic exponent a, but the position of the line 
has been chosen arbitrarily. 
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TABLES 



k 


7 


D2 


a 


Di 


^(1/2) 


^(3/4) 


1 


0.47 


0.49±0.01 


0.736±0.001 


0.7±0.11 


0.94 


0.91 


9 


0.53 


0.544±0.005 


0.755±0.001 


0.7±0.05 


0.96 


0.9 


16 


0.54 


0.549±0.004 


0.76±0.002 


0.75±0.05 


0.94 


0.89 


110.2 


0.564 


0.565±0.004 


0.76±0.002 


0.76±0.03 


0.94 


0.9 



TABLE I. Some Fractal Dimensions and Dynamical exponents . The ratio of the size of the 



interval / to 27r was the Golden Ratio for the 1st row, 2/3 for the others. 
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